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£> ' Abstract 

<n, 

\^ i We give large deviation upper bounds, and discuss lower bounds, for 

^^ . the Gibbs-KMS state of a system of quantum spins or an interacting 

^O | Fermi gas on the lattice. We cover general interactions and general 

~l ■ observables, both in the high temperature regime and in dimension 

f- «) ' one 

^! 1 Introduction and statement of the results 

The study of large deviations of macroscopic observables plays a fundamen 
tal role in classical statistical physics, for example in the study of the equiva- 
lence of ensembles and in hydro dynamical limits. The large deviation princi- 
ple for Gibbs measures in classical mechanics, initiated in the seminal papers 
ED HO] is now a classical subject, see e.g. jTHl EHl El [131 dni H21 ESI 1^1 E01 - 
It is maybe surprising that, in comparison, very little is known about large 
deviations in quantum statistical mechanics. To our knowledge, results on 
large deviations have been obtained only for the fluctuations of the par- 
ticle density in J2J> Ior ideal fermionic and bosonic quantum gases, and 
in |15| . for dilute fermionic and bosonic gases in (using cluster expansion 
techniques). As we were completing this paper a preprint [2H appeared 
where large deviation results for observables that depend only on one site 
are established in the high temperature regime (using again cluster expan- 
sions). In this work, however, we consider general observables. Also, from 
a technical point of view, we do not use cluster expansions, but a simple 
matrix inequality, combined with analyticity estimates on the dynamics and 
subadditivity arguments. 



Other large deviation results in quantum mechanical models can be 
found in |IJ |SJ |2B1 EHl ESJ- I n the larger context of probabilistic results 
for quantum lattice systems we want to mention |17U18ll251l26j . about the 
Central Limit Theorem and the algebra of normal fluctuations, and JH1E1) 
on Shannon-McMillan type of theorems. 

For simplicity we will consider quantum systems of spins or fermions 
on a lattice. The statistical mechanics of spins and fermions is naturally 
expressed in the formalism of C*-algebras, which we will use throughout the 
paper (refer to [HI EB] for spin systems and to the recent [2] for fermions) . 
For classical Gibbs systems, the DLR condition is a crucial ingredient in 
establishing the large deviations principle. For quantum systems, Araki 
introduced an analogous condition, the so-called Gibbs condition, which 
will play an important role in our proofs. Our results presumably extend 
to some bosonic systems or lattices of oscillators, but such systems present 
more technical difficulties as they are most naturally expressed in the W*- 
algebraic formalism. 

In order to illustrate the scope — and limitations — of this work, let us 
consider an example of a system of (fermionic) particles of spin 1/2. Our 
results, however, are more general and will be detailed in Section |3J 

Let c xa and c* a denote the annihilation and creation operators for 
fermions of spin a € {Til} a t site x € Z . These operators satisfy the 
Canonical Anticommutation Relations. We denote by n x%a = c* a c x%a the 
operator for the number of particles in x with spin a, so that n x = n x ,}-\-n x j, 
indicates the operator for the total number of particles in x. The finite- 
volume Hamiltonian (with free boundary conditions) for a finite subset A of 
the lattice is taken to be 

{x,y}cA cr X £A 

+ Y^ Jx- y n x n y . (1.1) 

{x,y}cA 

Special cases of this Hamiltonian are the Hubbard models and the t J-models, 
which are widely used in applications. We define the number operator for a 
finite subset A of the lattice by 

AA = ]>>-• ( L2 ) 

We work in the grand canonical ensemble, whose finite-volume Gibbs states 



are given by 

m \ trC • e ~' 9 ( fl ' A ~^ JV ' A )l 

^ A l j ~ t r(e-/ 3 (^A-^A)) ' ^ 

where (3 is the inverse temperature and is the chemical potential jjl. We 
denote by u/^'W the Gibbs states of the infinite systems, which, roughly 
speaking, are constructed as limit points of uj\ as A f Z d . Mathematically 
the Gibbs states of the infinite system are characterized, equivalently, by 
either the variational principle, the KMS condition, or the Gibbs condition. 
Let us consider a local (microscopic) observable A, i.e., a self-adjoint 
operator depending only on even products of creation and annihilation op- 
erators (this is a natural limitation for Fermi systems; see Section [2J). For 
example, say that A is an even polynomial of the creation and annihilation 
operators. Assume that A depends only on the sites of X C Z rf , and denote 
by v x A the translate of the operator A by x € Z d . For a cube A we define 
the macroscopic observable K\ by 

X+a;cA 

In this way, | A| — 1 i^A represents the average of A in A (|A| denotes the 
cardinality of the set A). If we denote by 1b the indicator function of a 
Borel set BcR, then, for a Gibbs state uj^^\ 

PA (B) =JM(l B (\A\- 1 K A )) (1.5) 



defines a probability measure on R. Namely, (|1.5|) is the probability that, in 
the state cu^ ,fJ- >, the observable |A|~ 1 ETa takes values in B. Large deviation 
theory studies the asymptotic behavior of the family of measures p\ on an 
exponential scale in |A|. This asymptotic behavior is expressed in terms of a 
rate function I(x), which is a lower continuous function with compact level 
sets. 

Let C be a closed set. We say that we have a large deviation upper 
bound for C if 

limsupT-j-rlogw^ Uc (IM^Ka)) < - inf I(x) . (1.6) 

A/Z d l A l V V " X&C 

Similarly, if O is open, we have a large deviation lower bound for O if 

liminf -Llogu/^ (lo (\Ar l K A )) > - inf I(x) . (1.7) 

A/"Z d |A| V w V 1 J J - xdO 



One says that {pa} satisfies the large deviation principle if we have upper 
and lower bound for all closed and open sets, respectively. 

In order to study the large deviations for p\ as A /" Z d (along sequences 
of cubes), one considers the corresponding logarithmic moment generating 
function, defined as 

f(a) = lim -Llogu^Ce ^). (1.8) 

A/"Z d |A| 

The Gartner-Ellis Theorem (see, e.g., ^J), shows that the existence of 
/(a) implies large deviation upper bounds with a rate function I(x) that is 
the Legendre transform of /(a). One obtains lower bounds if, in addition, 
the moment generating function is smooth, at least a C 1 . If the moment 
generating function is not smooth, one has a weaker result: in 1)1.70 . the 
infimum over O is replaced by the infimum over Oflfi, where E is the set 
of the so-called exposed points (see ^J for details). 

Our results apply both in one dimension and at high temperature. In 
both cases the parameters (3 and \i are such that there is a unique Gibbs- 
KMS state u&'A. 

Dimension one. Let us assume that the lattice is one-dimensional and 
that the interaction has finite range, i.e., there exists an R > such that 
T x _ y and J x -y vanish whenever \x — y\ > R. Our core result is that, for any 
macroscopic observable K\ and all values of j3 and /x, the moment generating 
function f{a) exists and is finite for all a S R. Furthermore f{a) is given 
by the formula 

f(a) = lim T-rrlog ; -r-— : rr^- — , (1.9) 

A^z d |A| tr(e~^(- ff A-/^VA)) v ' 

which involves only finite-dimensional objects. 

As recalled above, if I{x) is the Legendre transform of /(a), the Gartner- 
Ellis Theorem entails the large deviation upper bounds with / as the rate 
function. As for the lower bounds, it is tempting to conjecture that the 
function f(a) is smooth, in one dimension. We have not proved it so far. 

It is instructive, at this point, to compare our results on quantum systems 
with their classical analogs. In the classical case, using the DLR equations, 
one shows that a formula similar to Eq. (|1.9|) holds, with the trace replaced 
by the expectation with respect to the counting measure. In that case, 
one sees that /(a) is simply the translated pressure corresponding to the 
Hamiltonians (3H\ — (3[aN\ — aK\. Therefore, classically, the smoothness of 
/(a) follows immediately from the lack of phase transitions in one dimension, 



together with the identification of Gibbs states with functionals tangent to 
the pressure ^3 [""""]. In the quantum case, K\ does not commute with 
H\ — nN\, in general, so the thermodynamic interpretation of the moment 
generating function is not obvious. This is the main difference. 

Such interpretation is possible, however, when K\ commutes with H\ — 
fj,N\. In our example (jl.l|) . H\ does commute with N\ (this is not hard to 
verify, using the CAR; see Q2.1|) ). so that we can fully treat the physically 
important large deviations in the energy and in the density. More in detail, 
if the pressure function for our system is defined as 

P(P,ri = lim -Llogtr(V^-^A)A (L10) 

A/Z d |A| V / 

then Q1.6|) - Q1.7"i hold for the energy (K^ = H\), with I{x) being the Legen- 
dre transform of 

a -» P (p - a, J^») " P(J3,U) ■ (1.H) 

They also hold for the number of particles (K\ = N\), and in that case 
I(x) is the Legendre transform of 

a«P^+^-P(/3, M ). (1.12) 

High temperature. For arbitrary space dimension we assume that that 
the interaction is summable: J2 x ez d \T x \ + \Jx\ < oo. Our main result is that 
there exist two constants, /?o (which depends only on the Hamiltonians H\) 
and ao (which depends only on the observable K\), such that the function 
f(a) exists for \a\ < ao, \/3\ < /?o ; and arbitrary /tER. Furthermore, in the 
special case in which the macroscopic observable is a sum of terms depending 
only on one site, ao can be taken to be infinity. Again, the function /(a) is 
also given by (|1.9|h 

This yields large deviation upper bounds for closed sets which are con- 
tained in a neighborhood of the average K = liniA/z d lAI -1 ^;'^'^^^^)- At 
high temperature, one expects /(a) to be smooth, in fact analytic, and this 
can be proved using a cluster expansion [2*2] , 

For the case of commuting observables we show that the moment gen- 
erating function exists for any a, provided \(5\ < /3q. It is known that, for 
sufficiently high temperature and any value of the chemical potential, there 
is a unique Gibbs state (see Theorem 6.2.46 of """]). Using this, we obtain 
a full large deviation principle for the particle number (or density). As for 



the energy, we expect /(a) to have a singularity at some a 7^ 0; at any rate, 
we have upper bounds for all closed sets and lower bounds for sets that are 
contained in a neighborhood of the mean energy. For both the energy and 
the particle density, the rate functions are again the Legendre transforms of 

Once again, the precise statements for a general quantum lattice system will 
be presented — and proved — in Sectional 

2 Quantum lattice systems 

We consider a quantum mechanical system on the d-dimensional lattice Z rf , 
as seen, e.g., in JH2J E3 El EH] for spin systems, and in [H] for fermions. 

2.1 Observable algebras 

We first describe quantum spin systems. Let H be a finite-dimensional 
Hilbert space. One associates with each lattice site x 6 Z a Hilbert space 
7i x isomorphic to H and with each finite subset IcZ d the tensor product 
space TCx = <S> x ex %:• The local algebra of observables is given by Ox = 
S(TCx), the set of all bounded operators on TCx- If X C Y, there is a 
natural inclusion of Ox into Oy, and the algebras {Ox} form a partially 
ordered family of matrix algebras. The norm-completion of the union of 
the local algebras is a C*-algebra denoted by O which correspond to the 
physical observables of the system. In particular we have that [Ox , Oy] = 
whenever Xr\Y = %. A state u is a positive normalized linear functional on 
O, i.e., uj : O — ► C, uj(1) = 1 and oj(A) > 0, whenever A > 0. The group 
Z d acts as a *-automorphic group on O: For x € Z d , v x (Ox) = Ox+x- 
A state is called translation invariant if u> o v x = lu for all x £ Z d and we 
denote by fli the set of all translation invariant states. The action of v 
is asymptotically abelian: therefore Q/ is a Choquet simplex and one can 
decompose a state into ergodic components (see 33 ). 

The structure of the algebra of observables for fermionic lattices gases is 
a little more involved, due to the anticommutativity properties of creation 
and annihilation operators (see |31E]]). We construct it as follows. 

Let I be the finite set that is supposed to describe the spin states of a 
particle. For X a finite subset of Z rf , Tx is defined formally as the C*-algebra 
generated by the elements {c* a , c x ^} x& x,a&i together with the relations 



a'} = 0. (2.1) 



c* 


Cy,a' J 


~ $x,y8a,a' 1 


C* 
^x,cr : 


Cy,a' I 


~ \C-x,a •> Cy,< 



The above are referred to as CAR (Canonical Anticommutation Relations), 
c* a and c X)CT are called the annihilation and creation operators and are taken 
to be mutually adjoint by definition. It is easy to realize that, as a vector 
space, 

T X = span {c|\ j(T1 c^ ■ ■ ■ <fe >ffm } , (2.2) 

where the span is taken over all (finite) sequences {Ocj>0j>ttj)}7Li in X x 
Ix{-,*} that are strictly increasing w.r.t. a predetermined order. IfXcY, 
there is a natural inclusion J-x C Jy-, and we define the fermionic C*-algebra 
T to be the norm-completion of Uxcz d -^x- 

Elements of J- localized on disjoint parts of the lattice do not necessarily 
commute (they might either commute or anticommute) and so J- is not 
asymptotically abelian. We have to restrict the class of allowed observables 
to a smaller algebra. Let us denote by the automorphism of J- determined 
by 0(c| CT ) = — cjj.0.. The observable algebra of a fermionic lattice gas O is 
defined to be the even part of J 7 , i.e., 

O = {AeT\Q{A) =A}. (2.3) 

Clearly, Ox = O (1 Tx is given by the same r.h.s. of (|2.2|) . restricted to 
m even. Hence [Ox , 3~y\ = whenever X C)Y = 0, which is the com- 
mutativity property we need. The algebra O is thus quasilocal and similar 
considerations as for quantum spins systems apply. 

Example 2.1 For quantum spin systems with spin 1/2, the Hilbert spaces 
7i x , x G Z d is isomorphic to C 2 . 

Example 2.2 For fermionic systems of particles with spin 1/2, for each x, 
the algebra generated by c* a and c XjCT is isomorphic to £>(C 4 ). 

2.2 Interactions and macroscopic observables 

An interaction $ = {4>x} is a map from the finite subsets X of Z rf (de- 
noted Vf(Z d )) into the self-adjoint elements of the observable algebras Ox 
(denoted O^ ). We will always assume the interaction to be translation 
invariant, i.e., v x (j)x = <j>x+x f° r all x £ Z and all AT G Vf(Z ). An inter- 
action is said to have finite range if there exists an i? > such that </>x = 
whenever diam(X), the diameter of X, exceeds R. (One usually says that 
the range is R if R is the smallest positive number that verifies the previous 
condition.) We denote by B"' the set of all finite range interactions. The 



set of interactions can be made into a Banach space by completing B"' with 
respect to various norms. In this paper we use the norm 

MX = £ ||0x|| e A l x l, (2.4) 

XBO 

where A > and |X| denotes the cardinality of X. We call B\ the corre- 
sponding Banach space of interactions. To a given $ one associates a family 
of Hamiltonians (or energy operators) {H\} AeV , Z d\ via 

H A = ff A ($) = ]T <f>x- (2.5) 

XcA 

As in j20j , we define a finite-range macroscopic observable K of range R 
to be a mapping K : V f (Z d ) — ► C (sa) such that 

1. K A+X = v x K A for all x £ Z d and for all A G V f (Z d ). 

2. i^AuA' = -Ka + ^A' if A. and A are at distance greater than R. 

The kind of example that we have in mind, and that covers most applica- 
tions, is K\ = J2x+xgA v xA, for a given self-adjoint A £ Ox (which could 
be, say, the magnetization or the occupation operator at the origin, or the 
energy in a finite region, or so). 

Given a finite-range observable K, we can recursively define a finite-range 
interaction \I> G B^> by means of the equalities K\ = X^xcA ^x ■ We have a 
one-to-one correspondence between finite-range macroscopic observables and 
finite range interactions. We can and will consider more general macroscopic 
observables by replacing condition 2 with the condition that the interaction 
^, corresponding to K, belongs to some Banach space. 

2.3 Gibbs-KMS states 

There are several equivalent ways to characterize the equilibrium states cor- 
responding to an interaction <1>. These equivalences certainly hold if <£ G B\, 
for some A > |33l 15] . A more general result of this type has been proved 
recently in [3], both for spin and fermion systems, for a nearly optimal class 
of interactions. 

In this paper, the notation A / Z will always mean that we take the 
limit along an increasing sequence of hypercubes A. All our results can 
presumably also be proved for more general sequences (Van-Hove limits), 
but, for simplicity, we will refrain from doing so. 



We denote by -P(3?) the pressure for the interaction <£, given by the limit 
-P(3>) = lim A y. Z d |A| -1 tr(e~^ A ). Here tr is the normalized trace in ri\ 
and H\ is specified by (|2.5[) . Let w be a translation invariant state. The 
mean energy relative to u is defined as e™(w) = H = lim A * Z d \A\~ 1 lj(H\). 
Denoting by uj\ the restriction of oj to 0\, we define the mean entropy in 
the state u> by s(u>) = lim A ^ Z d \A\~ 1 S(u>a), where S(uja) = w^logp^) = 
tr(pA log pa) and p\ is the density matrix of uj\. The existence of the limits 
for the pressure, mean energy and entropy is a standard result. 

The variational principle states that 

P($) = sup (s{u) - e^\to)) . (2.6) 

We denote by Qj the set of states for which the supremum in Eq. (|2.6|) is 
attained, and we call such states the equilibrium states for the interaction 3>. 
The set (!j isa simplex and each of its states has a unique decomposition 
into ergodic states. 

The second characterization of equilibrium states is via the KMS condi- 
tion. Let us consider Tt, a strongly continuous unitary action of R on O. 
It is known that, on a norm-dense subalgebra of O, Tt can be extend to a 
(pointwise analytic) action of C (Hj. So, a state uj is said to be r-KMS if 

u(An{B)) = lo(BA) (2.7) 

for all A, B in a norm-dense r-invariant subalgebra of O. For a given 
interaction <J>, one constructs the dynamics r t as the limit of finite volume 
dynamics defined, on a local observable A, by e lHAt Ae~ lHAt . Then one can 
speak of a KMS state for the interaction <E>. 

The third characterization is through the Gibbs condition. This con- 
dition is analog to the DLR equations for classical spin systems. Stating 
it properly would require considerable machinery, including the Tomita- 
Takesaki theory. Detailed expositions can be found in [SJ ESI an d we will 
be brief here. Given an element P € 0^ sa > and a state to, one can define a 
perturbed state lo p in the following way: Using the Tomita-Takesaki the- 
ory one constructs (in the GNS representation) a dynamics Tt that makes 
u a r-KMS state. One then perturbs the dynamics Tt by formally adding 
the term i[P , •] to its generator (this would correspond to adding P to the 
Hamiltonian). Finally, one defines uj p as the KMS state for the perturbed 
dynamics (Araki's perturbation theory). 



For an interaction <3?, let us consider the perturbation 



W A = Y, <t>*> (2-8) 

which is well-defined under our assumptions. The state u) satisfies the Gibbs 
condition if, for every finite subset A, there exists a state u/ on C?a c such 
that 

co-Wa = J*) g, j _ (2.9) 

Here 0\c is the subalgebra of observables that "do not depend on A" (we 
omit the formal definition; suffices to say that O = 0\ <S> Oa c )- Also, which 
is crucial, oj a is the finite-volume Gibbs state on 0\ given by 



tr(e" 



"PiA) = ^ ' , (2.10) 



The Gibbs condition is very similar to the DLR equations in classical lattice 
systems, and it is not difficult to check that the DLR equations and the 
Gibbs condition are indeed equivalent for classical spin systems. 

Nor is it hard to verify that finite-volume Gibbs states satisfy all the 
previous three conditions. A fundamental result of quantum statistical me- 
chanics, due to Lanford, Robinson, Ruelle and Araki, asserts that the three 
characterizations are indeed equivalent for infinite- volume translation invari- 
ant states of spins or fermions. The key to the proof is the Gibbs condition, 
introduced by Araki. In the very recent [5|, equivalence has been proved 
for a very large class of interactions, much larger than the one considered in 
this paper. 

3 Moment generating function 

Given an interaction $ with a corresponding Gibbs-KMS state w G fil , 
and a macroscopic observable {i^A}> uniquely determined by the interaction 
\&, we introduce the moment generating function 

/(*'*>(£*) = lim -L\ogu>(e aK *); (3.1) 

A/Z 1 * |A| 



that is, when the limit exists. A priori it is not obvious that /( ' '(a) 
depends only on <3? and not the choice of lo G S7j . In this paper, however, 
we will always work in the one-phase regime, see Remark EH Furthermore 

10 



one expects that, as in the classical case, f^^'(a) would depend only on 

We will make one of the following assumptions. 
HI: High temperature. Both $ and ty belong to some B\ and 

J|%<1. (3-2) 

H2: High temperature improved. <3? is the sum of two interactions, <£ = 
<&' + $", where <J>" = {<p'x} x ez d involves only observables depending on 
one site, and, for all A C Z d , we have [H' A , H'^\ = 0. Also, we assume 
that $' and *$> belong to some B\ with 

^II*'IIa<1. (3-3) 

No smallness assumption on $" is made. 

H3: Dimension one. The lattice has dimension one and both $ and $> 
have finite range R. 

Remark 3.1 Condition H2 is important in physical applications where <£" 
is a chemical potential or an external magnetic field. It allows us to prove our 
results at high temperature for any value of the chemical potential/magnetic 
field (see the example in the introduction). 

Our main result is 

Theorem 3.2 Let u, <3? ; ^> be as above. 

1. (High temperature) // HI or H2 is satisfied, then the moment 
generating function /(*'*) (a) exists and is finite for all real a such 
that 

M<xmr <3 - 4) 

If the macroscopic observable is the sum of observables depending only 
on one site, i.e., K\ = J^xeA^xt withtp x € 0{ x y, then /(*>*) (a) exists 
and is finite for all a € R. 

2. (Dimension one) 7/H3 is satisfied, then /(*'*) (a) exists and is finite 
for all a € R. 



11 



The moment generating function /(*'*) (a) is convex and Lipschitz con- 
tinuous; more precisely, 

/(*.*>( ai ) - /<*'*> (a 2 ) | < ||*||oK - a 2 \ , (3.5) 

where ||\I/||o = J2x^o IIV'.X'IIJ- Moreover 

/ <^ (a)= lim _i log f!^^iy ( 3.o) 

J v ; A/z<* |A| ^ tr(e~ H A) y v ; 

Remark 3.3 Although our proof does not directly use this fact, the as- 
sumptions of Theorem 13.21 imply that there is a unique KMS state (in [S] , 
for instance, check Theorem 6.2.45 for HI, Theorem 6.2.46 for H2, and 
Theorem 6.2.47 for H3). 

Remark 3.4 The equality of the two limits (|3,1|) and H3.6J) implies that — 
using the terminology of [2j| — semi-local large deviations are the same as 
global large deviations. In other words, u;(1#(|A|~ 1 A'a)) decreases at the 
same exponential rate as wa(1b(|A|~ 1 A'a)). Global large deviations are so 
named because they gauge the probability of deviation from the expected 
value when a microscopic observable is averaged over all the available vol- 
ume. 

For particular, physically important observables, the results of Theorem 
2 can be improved. 



Corollary 3.5 Suppose that, for all A G Vf(Z d ), the observable K A com- 
mutes with the energy H\. 

1. If HI or H2 holds, then p ■ ■ ' >(a) exists and is finite for all a £ R, 
and is C 1 in a neighborhood of 0. // K\ is the sum of observables 
depending only on one site, then /' • \a) is C 1 for all a. 

2. If US holds, then f( ' '(a) exists, is finite, and is C 1 for all a E R. 
Proof: If [H A , K A ] = then, by Theorem E21 and Eq. (l3~6l) . 

/<«>(«, = lim 1 Ar (e ^-«-)\ 

J V ; A/zd |A| B \ tr(e-^A) J 

= P($ - a¥) - i>($) , (3.7) 
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so that, as in the classical case, /' ' '(a) is the translated pressure. There is 
a unique Gibbs-KMS state for the interaction <3? — a\&, provided ||<3? — c^Ha 
is sufficiently small (jH], Theorem 6.2.45), so, by the equivalence between 
Gibbs-KMS states and functionals tangent to the pressure HOES], /(*'*) (a) 
is differentiable if a is sufficiently small. If the interaction \& consists only of 
observables depending on one site, and H\ commutes with K\, then there is 
a unique Gibbs-KMS state for $ — a^f, for all a, provided \\&\\\ is small (jH], 
Theorem 6.2.46). If condition H2 is satisfied, similar considerations apply 
(see jH], Theorem 6.2.46). If condition H3 is satisfied there is a unique 
Gibbs-KMS state for $ - a® ( 8_, Theorem 6.2.47). I 

The proof of Theorem 13.21 is in two steps. In the first step, instead of 
/(*'*) (a), we consider 

5 (*. $ )(a) = lim -Llogtr(e aA ' A e- HA ). (3.8) 

A/*Z d |A| 

In the second step we show that 

/(*.*)(«) = 5 <*'*)(a)-P($). (3.9) 

The function g^^'{a) is defined via finite-dimensional objects. We will 
prove the existence of the limit using a subaddivity argument, as in the 
proof of the existence of the pressure. The equality ()3.9|) is proved using 
perturbation theory for KMS states. 

3.1 Perturbation of KMS states 

A basic ingredient in the proof of the existence of the pressure is the following 
matrix inequality: 

H+P\ i_ + „/0,H 



logtr(y +/ J -logtrfe^J < ||P|| , (3.10) 

where H and P are symmetric nxn matrices. In order to study the function 
g(*>*)( a ) ; where we have two (generally non-commuting) exponentials under 
the trace, one needs to estimate quantities like 



logtr(Ce^ +p ) -logtr(Ce H ) , (3.11) 



where C is a positive-definite nxn matrix. A little thinking convinces one 
that an estimate of (|3.11|) by a constant times ||P|| cannot possibly hold 
true, if the constant is required not to depend on C or n. 

The following lemma gives an upper bound for ()3.11|) which is indepen- 
dent of C and n, although it has a different form than Eq. H3.1U|) . 

13 



Lemma 3.6 Let H,P £ C nxn , with H* = H and P* = P. 

1. We have 

logtr (e p+p ) - logtr (e^) 

2. Also, ifCe C nxn with C > 0, 



< IIPII 



(3.12) 



logtr (Ce H+p ) - logtr (Ce H 



where 



< sup sup \\U- s {t)PU s {t)\\ , 

0<t<l _1< 5 <1 

(3.13) 

U s (t) = e< H+tp ^> . (3.14) 

Proof: The proof of part 1 is standard. One writes 



logtr (e H+p ) - logtr (e H 



1 dt4- logtr (c" +ir 
o at 



< / ds 
Jo 



tr(Pe 



H+tP\ 



ti(e H+tp ) 
having used the fact that, for E > 0, 

tr(AE) 



< \\P\ 



tr(P) 



< \\A\ 



(3.15) 



(3.16) 



To prove part 2, we recall DuHamel's identity for the derivative of e '", 
when F(t) is a bounded operator: 



d 

—i 
dt 



F « = / due uF ^F'(t)e^-^ F ^ 



We write 



logtr (Ce p+P ) -logtr (Ce^) = f dt — logtr (Ce 



H+tP 



(3.17) 



(3.18) 



and 



-logtr(Ce H+ < p 
dt V 



ti ^duCe u{H+tp )Pe 



(l-u)(H+tP) 



< 



ti(Ce H + tp ) 

tr L(H+tP)/2 C e (H+tP)/2 jl du e (u-l/2)(H+tP) p e (l/2-u){H+tP) 
tT ( e (H+tP)/2c e (H+tP)/2) 

[ 1/2 dse-< H+tp ^Ce^ H+tp ^ , (3.19) 

J -1/2 
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where we have used the bound (J3TR)Jl with E = e ( - H+tp ^ 2 Ce^ H+tp ^ 2 . This 
concludes the proof of Lemma 13, 61 | 

Lemma 13.61 involves the quantity U~ s (t)PU s (t), which is the time evo- 
lution (in imaginary time) of the observable P, relative to the dynamics 
generated by H + tP. One needs to estimate the dynamics for imaginary 
times between —i/2 and i/2. The connection with the KMS boundary con- 
ditions is evident. 

If we define a (finite-volume) state u) and a perturbed state u by 

,,, ti(Ae H ) PfA , tr(Ae H+p ) . . 

then Lemma 13.61 immediately implies that, for C > 0, 

log uj p {C)- log uj(C) < \\P\\+ sup sup || U- s (t)PU s (t) || . (3.21) 

0<t<l _I< S <I 

We will generalize this bound for Gibbs-KMS states of the infinite system, 
using results from the perturbation theory of KMS states (see, e.g., Chapter 
5.4 of .8; or Chapter IV. 5 of [53]). For a r-KMS state u, we denote by 
{Quj-, t^ui, OJ) = (G, ir, O) its GNS representation. The scalar product on Q is 
indicated with ( • , • ). For any i^Owe have 

w(A) = (0,ir{A)0) (3.22) 

and the dynamics r is implemented by some self-adjoint operator H on Q: 

vr(r s (,4)) = e isH k{A) e~ isH . (3.23) 

From now on we will identify an element A with its representative it (A). 
This is possible since the two-sided ideal {A € 0\uj{A* A) = 0} is trivial 
(|55j. Theorem IV. 4. 10), therefore it is the left multiplication on O (33 , 
Theorem 1.7.5). 

For P e 0^ sa \ t p {A) given by 



r p (A) = r s (A) 






00 PS PS\ 

+ Vj" dsi ds 2 - 

n=l J ° J ° 


/•Sn-1 

••/ ds n [r Sn (P),[-. 
Jo 


■IT S1 (P),T S (A)}]---}- 

(3.24) 



defines a strongly continuous semigroup of automorphims of O implemented 
by H + P: 

r p (A) = e »(H+p) Ae -i.(H+p) {325) 
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Moreover we have 

if (A) = if r s (A)(Tf )* = if r 8 (A)(Tf y 1 , (3.26) 

where the unitary operator 

if = e is{H+P) e -isH (3>27) 

has the following representation as norm-convergent series: 



vjw . s , S1 rSn-l 

rf = l + J2i n dsx ds 2 --- ds n T Sn {P)---r Sl (P) 

~L JO JO JO 



(3.28) 



Furthermore, f(s, P) = if O, defined on R, extends to a holomorphic func- 
tion f(z, P) on {z € C | < Imz < 1/2} (i.e., the function is continuous and 

bounded on the close strip, and analytic on its interior). In particular, O 

^p 

j/2' 



belongs to the (maximal) domain of if 2 , so that one can set 



O p = Tf /2 = e -(^)/ 2 e H l 2 O . (3.29) 

Araki's perturbation theory asserts that the state ui given by 



s<A) - l ° p ' A ° P) - (°- (r -v M(r ^°; , , ,„, 

() " (° P '° P > " (o,(rf /2 )-(rf /2 )o) {3W) 

is a r p -KMS state. 

The bound in Lemma 13.61 involves the norm of the imaginary-time evo- 
lution of the perturbation P. Therefore, for infinite systems, we will assume 
that P is an analytic element for the dynamics in the strip {|Imz| < 1/2}: 
by this mean that r z (P) extends to a holomorphic function in the strip, in 
the sense specified above. This is clearly a strong assumption and the main 
limitation of our approach. 

Theorem 3.7 Let uj be a t-KMS state and let P € O be a self-adjoint 
analytic element in the strip {\Im.z\ < 1/2}. Then, for all positive C € O 
we have 

logw p (C)-log W (C) < ||P||+ sup sup \\t£{P)\\ (3.31) 

o<t<i _i< s <! 
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Proof: The proof of Theorem 13,71 follows closely the proof of Lemma 13.61 
We first assume that C 1 ' 2 is an analytic element for the dynamics r — such 
elements form a dense subalgebra of O ( |33). Proposition IV. 4. 6). Rewriting 
Eq. (EHH1) as 

Tf = 1 + V (is)" / dm / du 2 • • • / du n r SUn (P) • • • r SUl (P) (3.32) 
~L JO JO Jo 

and recalling the hypothesis on P, it is easy to extend Tf to a holomorphic 
function on {|Ims[ < 1/2}. In light of Eq. (|3.26|) . then, we conclude that 
C 1 2 is an analytic element for r in that same strip, for all < t < 1. 
Using Eq. (J3.3UJ) we have 

\ogu p {C)-\oguj{C) = / olt — log iu tp (C) 

Jo at 

= jf ' * | [log (o , (r*f 2 )* cr^f 2 o) - log (o^ , o^)] . (3.33) 



We now claim that 






^ d,rf(P)r^ 2 . (3.34) 



Verifying (|3.34j) would amount to a simple application of DuHamel's formula 
()3.17|) . if H were a bounded operator. In the case at hand we need to work 
a little harder, even though we use the same idea. For e > 0, let II e be 
the projection on the invariant space of H defined by values of its spectral 
measure in [— e, e]. Then TT e = 1 — Il e is the projection on the orthogonal 
space. Set 

H e = U e HU e , H' e = TT e HTl! el P e = n e PII e . (3.35) 

Clearly, H e and P e are bounded operators and [H' e , H e ] = [H' e , P e ] = 0. By 
means of IJ3.17JI . and after a change of variable, we verify that 



± e -(H e +tP e )/2 e H e /2 



/ dse~ s(He+tPe) P e e s(He+tPe) e~ {He+tPe)/2 e He/2 . 
dt~ Jo 

(3.36) 
Now we multiply each factor above by the corresponding term e uHe (u = 
±1/2, ±s); these terms commute with everything. We obtain 

A fl/2 

| r V2=-/ dsr^{P e )Tf 2 . (3.37) 

17 



That (|3.37|) becomes 1)3.34(1 . as e — > +00, follows from (|3.24() — or rather 
its analytic continuation — and (|3,32l) . since P e is entire analytic for r s , and 
||P e -P|| -+0. 

Once 1)3.34)1 is settled, we can write 



d_ 
dl 



log <o,(r&)*cr&o 



u> 



tP frO 



■1/ 



a cfa7^(p)c+/^cfaC7ir(i') 



1/2 



The symmetric form of the KMS condition for 10 is easily derived from 
1)2.7)1 : for A, B analytic in the strip, 

J p (r%(A) r$(B)) = J p (BA) . (3.39) 

Applying the above twice, 

u> tP (ri P (P)C) = J P {r%{C^)ri( s+1/2) {P)r^C^) 



.tP 



We thus turn (l3~3"%l) into 



^' l r -j/2l° ) T Hs-l/2)K jr J''i/2 



(3.40) 



r/ 



itP 



s log(0-, CO 



tP 



u tp [r tp 



-i/2\ 



{C 1 ' 2 ) I%dsri P (P)rl P {C^) 



and therefore 



^lo K (O tp ,CO tp 
at 



utP (rtP /2{C i/i )T tP { ci/ 



< sup 



■5<«<5 



Here we have used the fact that 



A 



oj(B*AB) 



(3.41) 



(3.42) 



(3.43) 



uj{B*B) 

defines a state on O if uj(B*B) / 0. 

As for the second term in (J3.33)) . we plug C = 1 in 1)3.41)1 . use the 
invariance of a/ p with respect to r tp ', and conclude that 

d log(o tp ,O tFA 



dt 



< \\P\\. 



(3.44) 



This gives the desired bound when C 1 ' 2 is analytic. The general statement 
follows by density, see Corollary IV. 4. 4 in [33]. | 
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3.2 Analyticity estimates 

As is apparent from the previous section, we need estimates on the evolution 
of observables (in imaginary time). We will use two results, one valid at high 
temperature and one valid in dimension one. 

The first is due to Ruelle, has no restriction on the dimension, and is a 
standard. 

Proposition 3.8 Let $ G B\, for some A > (see Section WT^i . For any 
A G Vf(Z d ) and any collection of numbers {ux}xcA, with ux = "Ux(A) G 
[0,1], set 

H A ] = T. u x4>x (3.45) 

XcA 

(of course, H^ = H\, if ux = 1 for all X). If A G \J X Ox is a local 
observable and z belongs to the strip {|Imz| < 2/(A||$||a)}, then 



izH^ Ae -izH^ 



^ i , mi „ PI|e A|X| - (3-46) 



1 - \Jxnz\ % ||$||a 



This estimate is uniform in A (and {ux}) and thus holds in the limit 
A /" Z d , when this limit exists. In particular it holds for the infinite-volume 
dynamics t z . 

Proof: Follows trivially from the estimates of Theorem 6.2.4 in |S]. | 
Theorem 13.81 implies that, in the high-temperature regime 

^I|3|U<1, (3.47) 

local observables are analytic elements for the dynamics at least in the strip 
{|Imz| < 1/2}, which is what we need. 

The second estimate is due to Araki £Q and applies only in dimension 
one. It was used recently in [25^ to prove a central limit theorem in one- 
dimensional spin systems. 

In order to state it we introduce the concept of exponentially localized 
observables. Denote O n = 0[- n ,nl- Given A G O, we set ||-A||^ = ||-A|| and 

|| A|| W = inf \\A-A n \\. (3.48) 

This allows us to define, for < 6 < 1, the norm 

\\A\\ (g) = ^0- n |LA||M. (3.49) 

n>0 
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An element A of O is said to be exponentially localized with rate 9 if, and 
only if, 1 1 A 1 1(0) < oo. The symbol 0^> will denote the space of all such 
observables. 

We consider an interaction $ of finite range R, and set 



S(<&) 



■>K 



£ Q diam(A) 



(3.50) 



Also, for s > 0, we define 

Ms) = 

We have 



exp 



(-R+ 1)3 + 2^ 



R e kR_ 1 



fc=l 



(3.51) 



Proposition 3.9 Lei $ € i3^), wii/i range .R. If 9 € (0, 1) and ft, > verify 
Q e 4hS(®) _ 0/ ^ -^ f/ien i/iere exists a constant M = M(R, 9, h) (independent 
of $) such that, for A G 0^> and |Imz| < h, 



zli 



a ,4e~ 



izH 



(u) 



(6>) 



< MF R (2S($))\\A\ 



(9) 



(3.52) 



r(«) 



ifere ^^ a; is defined as in \S.4^jj - This estimate is uniform in A (and {ux}) 
and thus holds in the limit A /* Z , when this limit exists. In particular it 
holds for the infinite-volume dynamics r z . 



Proof: Follows from the results of PP; see also 



We will use this result in the particular case in which the macroscopic 
observable {K\} has finite range. Hence notice that, if A is a local ob- 
servable, then A € He C • Furthermore, for every 9 € (0,1), there exists 
a constant D = D(9,R') such that \\A\\^ < D\\A\\, for all A G O x with 
diam(A) < R' . The reverse bound, \\A\\ < \\A\\fg\, is of course valid for every 
A E O. These considerations and Proposition 13.91 imply that, for any such 
A, there exists a constant G = G(R, R' , S($>)) such that, for |Imz| < 1/2, 



izHa 



Ae 



-izH\ 



< G\\A\\. 



(3.53) 



Once again, this is what we need to apply Lemma 13.61 and Theorem 13. 71 
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3.3 Subadditivity 

We now give sufficient conditions for the limit 

g( 9 '*\a) = lim -L\ og tv(e aKA e~ HA ) (3.54) 

A/Z 1 * |A| 

to exist. 

Theorem 3.10 The following holds true: 

1. (High temperature) If condition HI or H2 applies, then the func- 
tion g( ' '(a) defined by Eq. \3. 5$ exists and is finite for a real, with 

\a\ < 1 „ T „ ■ (3.55) 

A||*|U 

Furthermore, if^ = {ipx} x &z d withip x G @{ x }> (observables depending 
only on one site), then g( ' '(a) exists and is finite for all a € R. 

2. (Dimension one) If condition H3 applies, then g^^> (a) exists and 
is finite for all a € R. 

Proof: We start with item 1 under the condition HI. The proof combines 
Lemma 13.61 the analyticity estimates of Section 13.21 and a subaddivity ar- 
gument as in the proof of the existence of the pressure. Let A to be an 
hypercube of side length L. We choose a > and write L = na + b, 
with < b < a. We divide the L-cube into disjoint adjacent n a-cubes, 
Ai, A2,- • • ,A n d and a "rest" region Ao which contains L d — {na) d lattice 
points. We write 

n d n d 

H A = Y, H A J +H Ao + W, K A = J2 K A j +K Ao +U. (3.56) 

3=1 3=1 

where 

W = Y,'<l>x, U = J2'^x (3.57) 

X X 

and J2 indicates a sum over all X C A such that, for some j = 0,1,... , n , 
X n Aj / and X n A] j£ 0. We denote by 

5f #) («) = ^logtr(e^e-^) (3.58) 
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the function whose limit we are set to take. By the commutativity property 
of local observables and the translation invariance, 



logtrfe^^i^.-E^i^ 



log Y[ tr (e""Je ' ~" 



aK A . -PH A , 



0=1 



= (na)*g% v '(a). (3.59) 

Set now 

P = H Ao + W , Q = K Ao + U. (3.60) 

Using Eq. (|3.59|) , the triangle inequality, Lemma 13.61 and Proposition 13.81 
we are able to estimate 



5a (»)--nr-9 Al («) 



< JL sup sup lle s(H A ~tP) pe s(H A -tP) l 



|A| 



0<t<l _I< S <I 

2- -2 



+ -|- SUp SUp \\ e -MKA-tQ) aTe sa(K A -tQ)\ 

|A| 0<<<1 -i< s <i 



< 



+ 



ml E +E' 

1 



i-^II*IIaIA| 



'xll e 



A|X| 



O' 



i-Mzll*IU l A l 



m E +£' ll«e- 



A|X| 



(3.61) 



! I! HA I I yxcAo X 

We take the limit A /* Z d of the various parts of Ineq. ()3.61j) . First, 



mE iM* A|Jf| < ^E El 

11 xeA xeA X9x 



*>x\\e 



L d -(na) d l .^. l 

< ^- $ A 



(3.62) 



as L — > oo; similarly for EjfeAo HV'xH e A ' x '. Also, in the same limit, 



^X>,ii« 



x\x\ 



< ^E E ll**|| e 



< 



D< 



n 



\\x\ 



xnA c ^0 



E 

xnA^« 



^ e 



A|X| 
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IAi 



E 



5x|| e 



a|x| 



(3.63) 



Once again, a similar estimate holds for Ylx IIV'xH e ' '. In the remainder, 
for the sake of the notation, we rename Ai = A. From (j3.6H) - (|3.63j) we 
obtain 



hm sup g A (a) - g y A (a, 

A/"Z d 



< 



I J_ y 

1 _ A||<b||. |A| ^ 

J- 4 II^IIA I"! xnA#0 
xnA c #0 



?X e 



A|X| 



+ 



a 



l-|a|lll*IUlA 



h £ 



,A[X| 



(3.64) 



xnA^B 
xnA c #0 



It is now time to take the limit A /* Z . Denote by A' the cube of side 
length a — a 1 ' 2 and concentric to A. We have 



A E 



3x|| e 



A|X| 



xnA^B 
xnA c ^0 



< 



< 



IAI 



IAI 



E E 

xeA' xsx 

xnA c ^0 

■ E 



ix||e A ^ + 



|A| 



E E 

ieA\A' x ^ 



, . A |v, |A\ A'L , 
l^xll e A|x| + ' X ' ||$| 



?x e 



0. 



A|X| 



(3.65) 



diam(X)>a 1 /2 

as a — > oo. The same holds for the second term of (|3.64l) . Finally, then, 



limsupgj* ' (a) — lim inf g\ ' (a] 

A/Z d A/Z d 



0. 



(3.66) 



which proves the existence and finiteness of the limit g( ' '(a), in the high 
temperature regime. 

In the special case in which \t consists only of one-body interactions, we 
have 






(3.67) 
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i.e., U = and all the observables involved commute. Thus, in the first 
inequality of ()3.61j) . the second term simplifies to 

\\ e - saiKA - tKA o) a K Ao e saiKA - tKA o)\\ = \\aK A \\ < {L d - (na) d )\a\ ||tf || . 

(3.68) 
Proceeding as above, one proves the existence of <?(*'*) (a) for all a£R. 

If condition H2 holds instead of HI we have to modify the argument 
a little: using the same notation as above and because <J>" only involves 
one-site interactions, we have 



H' A = J2H , A j +H' Ao + W, H'l = J2H Aj +H'i , (3.69) 

We note that since [H' v , H' v ] = for all V G Vf(Z d ), then the decomposi- 
tion l)3.69j) implies that 

[Hy, W] =0. (3.70) 

In order to estimate 



logtr (c< 



-Ha 



logtr 



f d 



(3.71) 



for positive C, we proceed in two steps, using Lemma l3~B1 We have, using 
(l3~70l . 



logtr (Ce- flA ) - logtr (Ce - ^ - ^ 



< sup sup 

0<t<l _I< S <1 

2— -2 



sup sup 

0<t<l _I< S <I 

2- s -2 

sup sup 

0<t<l -±< S <L 

2— S — 2 



e -s(H A -tW) W l e s(H A -tW) 
e -s(H' A -tW) e -sH'; w , e ,H» e s(H' A -tW) 
e - s (H' A -tW) w i e s(H' A -tW>) 



(3.72) 



This term does not involve <&" anymore and is estimated as under condition 
HI. On the other hand, since H\ — W' = Y^j^Aj + H Ao is a sum °f 
commuting terms, we have 



logtr (Ce"^-""') -logtr Ce"^^ 



< sup sup 

0<t'_ 

< \\H 



0<t<l _I< S <1 
2- -2 



{E^Ar'SJ^.^E^Ai-^ 



'All ' 



(3.73) 
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which is estimated as in (|3.68|) . 

If one works under condition H3 the proof is similar, using estimate 
(|3.53j) . This concludes the proof of Theorem I3.1UI I 



Theorem 3.11 If any of the conditions HI, H2, or H3 hold, and to is a 
Gibbs-KMS state for $ ; then 



lim 

k/z d 



i_ log .(^)_i_ los 



tr(e aKA e~ HA ) 



-H A ^ 



0. 



(3.74) 



Proof: We will give two different proofs of Theorem 13. Ill The first uses the 
Gibbs condition and not, a priori, the fact that we are in a one-phase region. 
Defining W\ as in (J2.8J) . we apply the Gibbs condition (|2.9|) for uj to the 
observable e aKA : 



UJ -W A ^aK^ = J*) ^aK^ ^ (1) 



tr (e aKA e~ H A 



(3.75) 



t r ( e -^A) 
On the other hand, Theorem 13 . 71 asserts that 

log uj~ Wa (e aKA ) -loguj(e aKA ) <||W A ||+sup sup 
v ' v ' 0<i<l -i< s <i 

(3.76) 

By proposition 13.81 then, if <J> is in the the high temperature regime and 

1*1 < 1/2, 



r- tWA (W A ) 



1 

W\ 



-tw A 



(W A 



< 



1 



1 



E 



ill^lUlAI xnA ^ 

xnA c ^0 



'xll e 



\\x\ 



(3.77) 



which vanishes when A /" Z d , as we have checked in ()3.65|) . The same, of 
course, happens to |A| _1 ||Wa||. Putting together ()3.75() . (|3.76|) and the last 
two estimates proves the theorem in the case HI. 

If H2 applies, we have only two relations that have to do with the specific 
case at hand: the rest of the proof are algebraic manipulations for KMS 
states. The first relation is 



•8 — ' S ~ '.S U 's 



(the notation should be clear), and the second is 

rf'\w K ) = W K . 



(3.78) 
(3.79) 
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Eq. ((5775)1 comes from the fact that [H' v , H v ] = 0, for all finite sets V C Z d , 
and that r s is the limit of finite- volume dynamics. As concerns Eq. ((3.79J) . 
we define 

w A (v) = Y, te- ( 3 - 8 °) 



xcv 
xnA^a ; xnA c #B 



As in the proof of Theorem 13.101 [H v , W A (V)] = 0, so, taking again the 
limit V / Z d , and noting that ||Wa(V) - W A \\ -> 0, we derive ((577J)) . 

Now, using ((3.78)1 and (|3.79() in (|3.24|) . we get that, for the perturbed 
dynamics, 

r- tw HA) = r f )'" iW/ A (rf ")(A)) . (3.81) 

We plug ^4 = Wa in the above, exploit ((3 .79)) again, and take the analytic 
continuation of the result: for \s\ < 1/2, 

r- tW \W k ) = rW<- tw \W A ) (3.82) 

which is estimated as in case HI. 

One proceeds similarly when H3 holds. This concludes the first proof of 
Theorem 13.111 

The second proof is based on the fact that — as we have thoroughly 
recalled earlier — the Gibbs-KMS state is unique, under our assumptions. 
Therefore we can write uj as limit of finite-volume Gibbs states with free 
boundary conditions: 

w(A) = lim tr( ^ e „ V } , (3.83) 

for A G Ux ®x- Let us write H v = H A + H A c + W A (V), where W A (V) was 
defined in Eq. (l3~50() . If A € O a , with A C V, 

tr{Ae- Hv ) tr{Ae- HA ) tr^e"-^--^-^^)) tr{e- H ^- H ^) 



tr(e- H v) tr(e-^A) tr^e"^"^) tT ( e -H A -H A c-w A (V)} ' 

(3.84) 
because the trace factorizes, when evaluating the product of two observables 
with disjoint support. Now, via Lemma f3.6[ a couple of estimates of the type 
seen in Theorem 13. 1UI yield 



lim TTT 
\yz d |A| 

lim — - 
A/Z d |A| 



tl f e -H A -H A . 

log 



log 



tr ( e -H A -H A c-W A (V)} 
tT r^ e -H A -H A o-W A (V)\ 



0, 



tr(Ae~ HA - H ^ c 
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(3.85) 



uniformly in A € Ox, A > 0, and in V D A. Thanks to this uniformity, one 
obtains the assertion of Theorem 13.111 from 1)3.84(1 . | 

We conclude by proving what we have called our main result. 
Proof of Theorem \3.2t Combining Theorems I3.1UI and 13.111 we have that 

/(*.*) (a) = lim -?-\ogLo{e aK ^) 
v ' A/Zd |A| v ' 

= lim — logtr(e aAA e _HA )- lim — logtr(e~ HA ) 
A/Z d |A| A/Z d |A| 

= 5 (*-*)(q)-P($). (3.86) 

The existence of the pressure is of course a standard result not harder than 
Theorem l3~TUl 

The convexity of /(*'*) follows from the convexity of a i— ► logw(e A ), 
which is verified with a standard application of Holder's inequality, noting 
that uj(e a A ) = / du(x) e ax , for some Borel measure i/ (coming from the 
spectral measure of K\ in the GNS representation). 

To obtain the Lipschitz continuity, we apply Lemma f3.6l with H = c%2K\, 
P = (oc\ — 02)ATa, and C = e~ A . Since H and P commute, 

j— logtr(e QlA ' A e- HA ) - logtr ( e ^ K ^ e ~ HA 

< T-rr\\(ai- a 2 )K A \\ < \ai-a 2 \^2yjx\\, (3.87) 

which easily leads to (|3.5|) . | 



Acknowledgments We thank Jean-Pierre Eckmann, Henri Epstein, Gio- 
vanni Gallavotti, Vojkan Jaksic, Claude-Alain Pillet, Charles-Edouard Pfis- 
ter, and Lawrence E. Thomas for their comments, encouragement, and sug- 
gestions. This work has its origin in the joint work of one of us (M. L.) with 
Joel Lebowitz and Herbert Spohn, to whom we are particularly grateful. 
L. R.-B. acknowledges the support from NSF Grant DMS-0306540. M. L. re- 
ceived travel funds from NSF (Block Grant DMS-0306887), G.N.F.M. (Italy) 
and COFIN-MIUR (Italy). 



27 



References 

[1] Araki, H.: Gibbs states of a one dimensional quantum lattice. Comm. 
Math. Phys. 14, 120-157 (1969) 

[2] Araki, H.: On the equivalence of the KMS condition and the variational 
principal for quantum lattice systems. Comm. Math. Phys. 38, 1-10 
(1974) 

[3] Araki, H. and Moriya, H.: Equilibrium statistical mechanics of fermion 
lattice systems. Preprint, 2003 

[4] van den Berg, M., Lewis, J. T., and Pule, J. V.: Large deviations 
and the boson gas. In: Stochastic mechanics and stochastic processes, 
Lecture Notes in Math. 1325, Berlin: Springer 1988, pp. 24-39 

[5] van den Berg, M., Lewis, J. T., and Pule, J. V.: The large deviation 
principle and some models of an interacting boson gas. Comm. Math. 
Phys. 118, 61-85 (1988) 

[6] Bjelakovic, I., Kroger, T., Siegmund-Schultze, R., and Szkola A.: 
The Shannon-McMillan theorem for ergodic quantum lattice systems. 
Preprint, 2003 

[7] Bjelakovic, I., Kroger, T., Siegmund-Schultze, R., and Szkola A.: 
Chained typical subspaces: a quantum version of Breiman's Theorem. 
Preprint, 2003 

[8] Bratteli, O. and Robinson, D. W.: Operator algebras and quantum 
statistical mechanics. Vols. 1 and 2. Texts and Monographs in Physics. 
Berlin: Springer Berlin, 1981 

[9] Comets, F.: Grandes deviations pour des champs de Gibbs sur Z . C. 
R. Acad. Sci. Paris Ser. I Math. 303, 511-513 (1986) 

[10] Datta, N., Fernandez, R., Frohlich, J., and Rey-Bellet, L. Low- 
temperature phase diagrams of quantum lattice systems. II. Convergent 
perturbation expansions and stability in systems with infinite degener- 
acy. Helv. Phys. Acta 69, 752-820 (1996) 

[11] Dembo, A. and Zeitouni, O.: Large deviations techniques and appli- 
cations. Second edition. Applications of Mathematics, 38. New York: 
Springer, 1998 



28 



[12] Deuschel, J.-D., Stroock, D. W., and Zessin, H. Microcanonical distri- 
butions for lattice gases. Comm. Math. Phys. 139 83-101 (1991) 

[13] Ellis, R. S.: Entropy, large deviations and statistical mechanics. 
Grundlehren der mathematischen Wissenschaften 271. Berlin Heidel- 
berg New York: Springer, 1985 

[14] Follmer, H. and Orey, S.: Large deviations for the empirical field of a 
Gibbs measure. Ann. Probab. 16, 961-977 (1988) 

[15] Gallavotti, G., Lebowitz, J. L., and Mastropietro, V.: Large deviations 
in rarefied quantum gases. J. Statist. Phys. 108, 831-861 (2002) 

[16] Georgii, H.-O.: Large deviations and maximum entropy principle for 
interacting random fields on Z d . Ann. Probab. 21, 1845-1875 (1993) 

[17] Goderis, D., Verbeure, A., and Vets, P.: Noncommutative central lim- 
its. Prob. Th. Rel. Fields 82, 527-544 (1989) 

[18] Goderis, D., Verbeure, A., and Vets, P.: Dynamics of fluctuations for 
quantum lattice systems. Comm. Math. Phys. 128, 533-549 (1990). 

[19] Israel, R. B.: Convexity in the theory of lattice gases. Princeton series 
in physics. Princeton: Princeton University Press, 1979 

[20] Lanford III, O. E.: Entropy and equilibrium states in classical statistical 
mechanics. In Statistical mechanics and mathematical problems Lecture 
Notes in Physics 20, Berlin Heidelberg New York: Springer, 1973, pp. 
1-113 

[21] Lebowitz, J. L., Lenci, M., and Spohn, H.: Large deviations for ideal 
quantum systems. J. Math. Phys. 41, 1224-1243 (2000) 

[22] Lenci, M. and Rey-Bellet, L.: Large deviations in quantum lattice 
systems: Analiticity at high temperature. In preparation 

[23] Lewis, J. T., Pfister, C.-E., and Sullivan, W. G.: The equivalence of 
ensembles for lattice systems: some examples and a counterexample. J. 
Stat. Phys. 77, 397-419 (1994) 

[24] Lewis, J. T., Pfister, C.-E., and Sullivan, W. G.: Entropy, concentration 
of probability and conditional limit theorems. Markov Process. Related 
Fields 1, 319-386 (1995) 



29 



[25] Matsui, T.: Bosonic central limit theorem for the one-dimensional XY 
model. Rev. Math. Phys. 14, 675-700 (2002) 

[26] Matsui, T.: On the algebra of fluctuation in quantum spin chains. 
Preprint (2002) 

[27] Netocny. K. and Redig, F.: Large deviations for quantum spin systems. 
Preprint (2004) 

[28] Ohya, P. and Petz, D.: Quantum entropy and its use Texts and Mono- 
graphs in Physics. Berlin: Springer, 1993 

[29] Olla, S.: Large deviations for Gibbs random fields. Probab. Th. Rel. 
Fields 77 343-357 (1988). 

[30] Pfister, C.-E.: Thermodynamical aspects of classical lattice systems. 
In In and out of equilibrium (Mambucaba, 2000), Progr. Probab., 51, 
Birkhauser Boston, Boston, MA, 2002. pp. 393-472 

[31] Ruelle, D.: Correlation functionals. J. Math.Phys. 6 201-220 (1965) 

[32] Ruelle, D.: Statistical mechanics: rigorous results River Edge, NJ: 
World Scientific, 1999 

[33] Simon, B.: The statistical mechanics of lattice gases Vol. I. Princeton 
series in physics. Princeton: Princeton University Press, 1993 

[34] Raggio, G. A. and Werner, R. F.: Quantum statistical mechanics of 
general mean field systems. Helv. Phys. Acta 62, 980-1003 (1989) 

[35] Raggio, G. A. and Werner, R. F.: The Gibbs variational principle 
for inhomogeneous mean-field systems. Helv. Phys. Acta 64, 633-667 
(1991) 



30 



